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BACKGROUND OF THE INVENTION 
1. Field of the Invention 

[0003] The present invention relates to a method for modeling a semiconductor 

device and more particularly to a method of modeling the thermal and electrical 
characteristics of a semiconductor device which utilizes a semiphyscial device model 
coupled with an analytical thermal resistance model to self consistently solve for the 
channel temperature and internal charge/electrical field structure of the semiconductor 
device. 

2o Description of the Prior Art 

[0004] HEMT technology provides RF components that have unparalleled, high 
performance characteristics at high frequencies (microwave to millimeter wave) and 
high power levels. As such, HEMTs are known to be used in various RF applications. 
Unfortunately, high power level applications also require high levels of DC power 
dissipation that elevates the HEMT components to high temperature levels. Currently, 
there are two predominant methods for modeling the electrothermal and thermal 
characteristics of HEMT device; finite element thermal simulation and physical device 
simulation 

[0005] Finite element thermal simulations are used to simulate the lay out 
dependent thermal conducting characteristics of a semiconductor device. The 
simulation may be accomplished either by simulating the two dimensional cross section 
of the device lay out and then assuming semi-infinite thermal conditions in the 
orthogonal direction; so called quasi-three dimensional modeling; or fully simulating 
the three dimensional device layout. Because the full three dimensional approach 
requires much more computational power and sophisticated software, the quasi-three 
dimensional approach is known to be in more common use. A typical example of this 
approach is shown in Fig. 1 which depicts a finite element mesh for a HEMT device 
layout. The device layout is typical of all HEMT devices. 

[0006] Finite element thermal simulations are known to provide accurate 
estimations of a devices thermal conduction. However, the main draw back with this 
method is the inability to couple the calculated channel temperature of the device back 
to an electrical simulation in order adjust the electrical characteristics of the device. In 
addition, this approach assumes a heat source, typically DC power dissipated within the 
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intrinsic device which remains constant. However, in reality the DC power dissipation 
also changes with temperature. Thus, for an accurate electrothermal simulation, the so- 
called self heating effect must be taken in account. 

[0007] As mentioned above, the full electrothermal characteristics can be 
simulated from so-called physical device simulators. Physical device simulators are 
known to utilize comprehensive knowledge about material characteristics and basic 
device physics to simulate the physical operation within the structure of the device. 
Usually simulators are based upon finite element or Monte Carlo approaches- although 
the finite element approach is usually employed for electrothermal simulations that also 
incorporate thermal conduction. Because these tools use the physical structure to 
simulate performance, the correspondence between the simulated electrothermal 
performance and the devices physical characteristics are relatively strong. However, 
the ability of device simulators to accurately model real, measured high frequency 
electrical characteristics is relatively inaccurate. Usually such simulation tools are able 
to achieve fairly useable modeling of DC characteristics but cannot be used for accurate 
high frequency simulations. Thus, there is a need for an electrothermal semiconductor 
device model which provides relatively accurate results at high frequencies. 

SUMMARY OF THE INVENTION 
[0008] Briefly the present invention relates to a method for modeling 
semiconductor devices which utilizes a semiphysical device model coupled with an 
analytical thermal resistance model to self consistently solve for the channel 
temperature and internal charge/electric field structure of the semiconductor device. As 
such the method in accordance with the present invention can realistically simulate the 
response of electrical performance to temperature and vice versa of a semiconductor 
device. 
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DESCRIPTION OF THE DRAWTNCS 
[0009J These and other advantages of the present invention will be readily 
understood with reference to the following specification and attached drawings 
wherein: 

[0010] FIG. 1 illustrates aknown finite element mesh for aHEMT device layout. 

[001 1] FIG. 2 is a graphical illustration illustrating the measured semi-physical 

modeled DC-IV characteristics at an ambrent temperature of 200° C. 

[0012] FIG. 3 is similar to FIG. 2 but at an ambient temperature of 25°C. 

[0013] FIG. 4 is a graphical illustration illustrating the measured vs semi- 

phyically modeled DC-IV characteristics at an ambient temperature of 25°C for a self 

heating large periphery 8 fingered 600^m gate periphery device cell. 

[0014] FIG. 5 is schematic diagram of an exemplary small signal equivalent 

circuit model for a HEMT device. 

[0015] FIG. 6 is a sectional view of an exemplary HEMT illustrating the rough 
translation of the physical origins for each of the equivalent circuit elements illustrated 
in the small signal circuit model in FIG. 1. 

[0016] FIG. 7 is a cross-sectional view of a HEMT illustrating the regions in the 
HEMT which correspond to the various circuit elements in the small signal equivalent 
circuit model illustrated in FIG. 5. 

[0017] FIG. 8 is an example of a relatively accurate measured-to-model I-V 
characteristics using the semi-physical modeling method in accordance with one aspect 
of the present invention. 

[0018] FIG.9 is a elevational view illustrating an epi stack for an exemplary 
HEMT, 

[0019] FIG. 1 0 is a cross-sectional view of a HEMT and an exemplary epi stack. 
[0020] FIG. 11 isablownupdiagramofthecross-sectionalparameterspertaining 

to the T-gate geometry for the exemplary epi stack illustrated in FIG. 7. 
[0021] FIG. 12 is a diagram of an electric conductance model used in the semi- 
physical example. 

[0022] FIG. 13 is a Smith chart illustrating the measured vs modeled S- 

parameters S11,S12 and S22 simulated in accordance with the method in accordance 
with the present invention. 

[0023] FIG. 1 4 illustrates the measured vs modeled values for the S2 1 parameter. 
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[0024] FIG. 15 is similar to FIG. 14 but for the S12 S-parameter. 

[0025] FIG. 1 6 represents an exemplary S-parameter microscope in accordance 

with the present invention. 

[0026] FIG. 17 illustrates the internal and external regions of an exemplary 
HEMT device. 

[0027] FIG. 18 is similar to FIG. 16 but illustrates the approximate locations of 

the model elements in the HEMT FET device illustrated is FIG. 16. 
[0028] FIG. 19 is a schematic diagram of a common source FET equivalent 

circuit model. 

[0029] FIG. 20 is an illustration of specific application of the S-parameter 

microscope illustrated in FIG. 16. 

[0030] FIG. 21 is similar to FIG. 16 which demonstrates the inability of known 

systems to accurately predict the internal charge and electrical field structure of a 
semiconductor device. 

[0031] FIG. 22 is a plan view of a four-fingered, 200^im GaAs HEMT device. 
[0032] FIG. 23 is a graphical illustration illustrating the measured drain-to-source 
current I ds as a function of drain-to-source voltage Vds for the sample FET device 
illustrated in FIG. 22. 

[0033] FIG. 24 is a graphical illustration illustrating the drain-to-source current 

I ds and transconductance G m as a function of the gate-to-source voltage V gs of the 
sample FET device illustrated in FIG. 22. 

[0034] FIG. 25 is a Smith chart illustrating the measured Sll, S12 and S22 
parameters from frequencies of 0.05 to 40.0 GHz for the FET device illustrated in FIG. 
22. 

[0035] FIG. 26 is a graphical illustration of the magnitude as a function of angle 

for the S21 S-parameter for frequencies of 0.05 to 40 GHz for the exemplary FET 
illustrated in FIG. 22. 

[0036] FIG. 27 is a graphical illustration of a charge control map of the charge 
and electric field distribution in the on mesa source access region shown with R s as a 
function bias in accordance with the present invention. 

[0037] FIG. 28 is a graphical illustration of a charge control map of charge and 
electric field distribution in the on-mesa drain access region shown with R d as a 
function of bias in accordance with the present invention. 
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[0038] FIG. 29 is a graphical illustration of a charge control map for the non- 
quasi static majority carrier transport, shown with 'R i as a function of bias in accordance 
with the present invention. 

[0039] FIG. 30 is a graphical illustration of a charge control map for gate 
modulated charge and distribution under the gate, shown with Cgs and Cgt as function 
of bias in accordance with the present invention. 

[0040] FIG. 3 1 is a plan view of an exemplary ir-FET with two gate fingers. 
[0041] FIG. 32 is a plan view of a tu-FET with four gate fingers. 
[0042] FIG. 33 is an illustration of a tt-FET parasitic model in accordance with 
the present invention. 

[0043] FIG. 34 is an illustration of an off-mesa parasitic model for a ti-FET in 
accordance with the present invention. 

[0044] FIG. 35 is an illustration of an interconnect and boundary parasitic model 
in accordance with the present invention for the tt-FET with four gate fingers as 
illustrated in FIG. 32. 

[0045] FIG. 36 is an illustration of an inter-electrode parasitic model in 
accordance with the present invention. 

[0046] FIG. 37 is a schematic diagram of the inter-electrode parasitic model 
illustrated in FIG. 36. 

[0047] FIG. 3 8 is an illustration of an on-mesa parasitic model in accordance with 
the present invention. 

[0048] FIG. 39 is a schematic diagram of the on-mesa parasitic model illustrated 
in FIG. 41. 

[0049] FIG. 40 is an illustration of an intrinsic model in accordance with the 

present invention. 

[0050] FIG. 41 is a schematic diagram of the intrinsic model illustrated in FIG. 

40. 

[0051] FIG. 42A is an exemplary device layout of a tu-FET with four gate fingers. 

[0052] FIG. 42B is an equivalent circuit model for the ti-FET illustrated in FIG. 
42A. 

[0053] FIG. 43 is a single finger unit device cell intrinsic model in accordance 

with the present invention. 
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[0054] FIG. 44 is similar to FIG. 43 and illustrates the first level of embedding 

in accordance with the present invention. 

[0055] FIG. 45 is similar to FIG. 43 and illustrates the second level of embedding 
in accordance with the present invention. 

[0056] FIG. 46 is an equivalent circuit model of the tc-FET illustrated in FIG. 
42A in accordance with the present invention. 

[0057] FIG. 47 is similar to FIG. 45 and illustrates the third level of embedding 
in accordance with the present invention. 

[0058] FIG. 48 is similar to FIG. 45 and illustrates the fourth level of embedding 

in accordance with the present invention. 

[0059] FIG. 49 is similar to FIG. 45 and illustrates the fifth level of embedding 
in accordance with the present invention. 

[0060] FIG. 50A and 50B is a flow chart of a parameter extraction modeling 

algorithm that forms a part of the present invention. 

[0061] FIGS. 5 1 and 52 represent the error metric in accordance with the present 
invention. 

[0062] FIG 53 A is a Smith chart illustrating the measured versus the initial model 

solutions for the Sll, S12 and S22 S-parameters from frequencies from 0.05 to 40.0 

GHz, 

[0063] FIG. 53B is a graphical illustration of angle versus magnitude for the 

initially modeled S-parameter S21 from frequencies of 0.05 to 40 GHz. 
[0064] FIG. 54A is a Smith chart illustrating the measured versus simulated S- 
parameters Sll, S12 and S22 for frequencies 0.05 to 40 GHz for the first extraction 
optimization cycle. 

[0065] FIG. 54B is a graphical illustration of magnitude as a function of angle for 
the measure and first optimized model S-21 parameter for frequencies 0.05 to 40 GHz 
for the first optimization cycle. 

[0066] FIG. 55 A is a Smith chart illustrating the measure as a function of the 
final model solution for S-parameters S 1 1 , S 12 and S22 for frequencies 0.05 to 40 GHZ 
for the final solution. 

[0067] FIG. 55B is a graphical illustrations of the magnitude as a function of an 
angle for S-parameter S21 for the final model solution from frequency 0.05 to 40 GHz. 
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[0068] FIG. 56 is a graphical illustration of the semi-physically modeled vs 
measured small signal Gm. 

[0069] FIG. 57 is a graphical illustration of the semi-physically simulated bias 
dependence of the small-signal output conductance Rds. 

[0070] FIG. 58 is a graphical illustration of the semi-physically simulated bias 
dependence of the small signal gate-source and gate-drain capacitance Cgs and Cgd. 
[0071] FIG. 59 is a graphical illustration of the semi-physically simulated bias 
dependence of the small signal gate source charging resistance Ri. 
[0072] FIG. 60 is a graphical illustration of the semi-physically bias dependence 
of the small signal source and drain resistance Rs and Rd. 

[0073] FIG. 61 is a graphical illustration of the measured vs modeled bias 
dependent gain at 23.5 Ghz for a K-band MMIC amplifier. 

[0074] FIG 62A and 62B are graphical illustration of the extracted parameters 
from measured device I-V's for process control monitor testing. 
[0075] FIG. 63 is a graphical illustration of the measured vs semi-physically 
simulated process variation for Gmpk and Vspk. 

[0076] Fig. 64 is a graphical illustration of the measured vs semi-physically 

simulated process variation for Idpk and Gmpk. 

[0077] FIG. 65 is a graphical illustration of the measured vs semi-physically 
simulated process variation for Imax and Vpo. 

[0078] FIG. 66 is a graphical illustration of the measured/extracted vs semi- 
physically simulated process variation for the small signal equivalent model Rds and 

Gm. 

[0079] FIG. 67 is a graphical illustration of the measured/extracted vs semi- 
physically simulated process variation for the small signal equivalent model Cgs and 

Gm. 

[0080] FIG. 68 is a graphical illustration of the measured vs semi-physically 

simulated physical dependence for Imax as a function of physical gate length. 
[0081] Fig, 69 is a graphical illustration of the measured/extracted model vs 

semi-physically simulated physical dependence for Rds as a function of physical recess 
undercut width. 
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DETAILED DESCRIPTION 
[0082] The model in accordance with the present invention utilizes a semi 
physical device model that is coupled to an analytical thermal conductance model as a 
means of simulating the electrothermal performance characteristics of a semiconductor 
device. In particular, a semiphysical model for HEMT devices is demonstrated that is 
able to relatively accurately represent small signal, noise, non-linear and large signal 
characteristics. In order to incorporate temperature dependence and also layout 
dependence of thermal conductivity, the following procedure is utilized: 
[0083] 1 ) Derive a Semi-Physical device model that is able to replicate measured 
DC I-V characteristics and bias-dependent small-signal characteristics 
very accurately at room temperature. 
[0084] 2) Incorporate any known temperature dependence of material parameters. 
[0085] 3) Measure DC-IV's and bias-dependent small-signal S-parameters across 

the desired range of temperatures. 
[0086] 4) Extract small-signal equivalent circuit models for each of the S- 

parameter measurement vs temperature. S-parameter microscopy as 
discussed below may be used to develop charge-control map solutions. 
[0087] 5) Develop temperature coefficient expressions that apply to the intrinsic 
device model expressions of the semi-physical model. These temperature 
coefficients must adjust the predictions of the semi-physical device model 
to match the measured DC and small-signal data at each temperature. In 
particular: 

[0088] a) I-V's must be matched. 

[0089] b) If conventional, prior art small-signal model extraction is 

performed, for example, Minasian extraction algorithms, as 
discussed in detail in "Broad Band Determination of the FET 
Small-Signal Equivalent Circuit", by Berroth, et al., IEEE-MTT . 
Vol 38, No. 7, July 1990. relative changes in the C-V's 
(capacitance-VS-voltage) must be matched by the semi-physical 
model predictions. 

[0090] c) If S-parameter microscopy is performed, absolute changes in 

the C-V's must be matched to the semi-physical model 
[0091] 6) Implement an appropriate analytical thermal conduction model. 
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[0092] 7) Couple the semi-physical device model and the analytical thermal 

conduction model by: 
[0093] a) Substituting the "environmental temperature" that operates in 

all of the temperature-dependent terms and temperature 

coefficients with "channel temperature" 
[0094] b) Using the semi-physically modeled length of the saturated 

region, XSAT, as the length of the heat-generating region. 

[0095] As will be discussed below, the semiphysical HEMT model may be 
extended to perform full electrothermal device simulations. In accordance with the 
above, a semiphysical model mentioned in step 1 is discussed below. In step 2, 
reported thermal properties of the material systems used for the HEMT technology are 
incorporated into an appropriate material related expression of the semiphysical device 
model These properties may be obtained from "GaAs, AlAs and AlxGal-xAs 
Materials Parameters For Use and Research and Device Applications" by S. Adachi, 
J.Appl. Phys, vol. 58, No. 3, August 1985, hereby incorporated by reference. In step 
number 3, the DC I-V and bias dependent S -parameters are measured for a standard 
device lay out across several base plate temperatures; for example -25°C, 25°C, 125°C 
and 200°C. In step 4, S-parameter microscopy is employed for each set of temperature 
dependent data. In step 5, a temperature compensation co-efficient is developed from 
reported material thermal dependance relationships, for example, as shown below. 



Ambient Temperature Coefficient TCamb [] = 1.247*1^ + 1.51^-0.41^ + (T aryf) - 300)(-3.95-1 . 15*^^*0.0001 
Room Temperature Coefficient TCref [] = 1.247^^^ + 1.51^- 0.4111/ + (273- 300)(-3.95-1. 15*^^*0.0001 
Temperature Coefficient Factor TCF [] =TC ani5 /TC ref 



[0096] Once the expression, TCF is applied to the appropriate semiphysical 
charge and transport control equations, temperature dependent I-V and C-V 
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characteristics are accurately fit. As shown, the equations below illustrate how key 
charge control and carrier transport relationships are enhanced to incorporate 
temperature dependence through the empirical temperature compensation term set forth 
below: 



l^maxO ^ "maxL v ds 

Maximum Channel Charge [cm" 2 ] H chan / H chanREF ) * TC F 

=(N S ' + N s )*(dj + Ad^ MshKpar * V £ 
Parasitic Charge Accumulation in the Donor Layer N sDonor [cm -2 ] (d; + Adj - * TC F ) 

Saturation Velocity V s [cm/s] 1 ,47E+07*TC F 

[0097] FIG. 3 demonstrates the enhanced semiphysical device model now able 

to accurately simulate I-V characteristics at 200°C, while Figure 4 demonstrates the 
accuracy of the I-V simulation at 25°C. Thus, the semiphysical model, which can be 
derived at room temperature can be modified with empirical temperature coefficients 
to incorporate accurate temperature dependence. 

[0098] In step 6, the analytical thermal conduction expression is as set forth in 
"Precise Technique Finds FET Thermal Resistance" by H. Cooke, " Microwaves and 
RF", August 1986. Lastly, in step 7, gate length, which is taken to be the length of the 
heat-generating region in Cooke's expressions, is replaced by a semiphyically modeled 
XSAT expression provided below. 
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V mM+rt/ A/ \«ni(1/im*1>i l 


Position of the Boundary between Regions 1 and 2 














= 5X { 1 + [Lgeff-ALs-Xs-Lg+X,^]^ 5 X) 








+ sqrt ( 5 2 + ([L flef rAL s -X s -L g +X satREF ]/(2 5 X) - 


Length of the Saturated Region, Region 2, under the Gate 


^SAT 


tun] 


} 


Length of the Linear Region, Region 1 




turn] 


= l-geff" ^SAT 








= S X{1+[L 98 -L 1 y(2 5 X) 


Total Length of the Saturated Region, Region 2 


L 2 


[n"i] 


+ sqrt( 5 2 + ([L gs -L 1 ]/(2 5X)-1) 2 } 



[0099] In Cooke's original derivation, it is assumed that the heat-source within 

the FET can be approximated by a uniform source for the length of the physical gate 
length. In reality, a more accurate expression for this heat source would have a length 
equal to the saturation region's length, or XSAT. Most of the drain voltage, and hence 
most of the DC power dissipation is dropped across this region, thus making this 
dimension more suitable to describe the source of heat. 

[01 00] The temperature compensation coefficient is modified to operate between 

the "channel temperature" rather the "ambient temperature" as generally shown by the 
equations provided below. 



Artiei ToTinduBCMidei: TCfchsn [] 
RrmTarpa^uBCMidal TChef [] 
Tet|H^ueCMicia1Fa±r TCF [] 
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[01011 The distinction between TCF and TCF is that the TCF is a temperature 
coefficient determined through complete, uniform heating of the device sample. When 
"base-plate" heating is used to heat the device sample, it can be assumed that base-plate 
temperatures are roughly equal to the "channel temperature", or rather the temperature 
of the device is the same as the surrounding environment. 

[0102] The implementation of Cooke's method coupled with the semi-physical 

device model results in the equations below: 
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Channel Temperature T CH [°CJ = e P CHSAT W g /1000 + T sink - 273 

= N f / { 27c[(N r 1)/ln(M)-(N r 2)/In{P)} [ K TH W g /1e7 
Thermal Conductivity © [°C/mW] ] } 

Power Dissipation over the Saturated Region P qhsat [mW/mmS = l dsW * V CHSAX * 10 6 



Cooke's "M" parameter 

Cooke's "P" parameter 
Thermal Resistance 
Drain Current / unit Gate Width 



M 

P 
KTH 

' dsW 



[] 

[] 

[W/cm°CJ 



= { 2 sqrt { cosh [n(d ff 
cosh W ffd -X SAT )/(4H sub )]} + 1}/ 
{sqrt { cosh [«(d ffd + X SAT )/(4H sub )] 
cosh [^(d ffd -X SAT )/(4H sub )]}-1 } 

2 sqrt {[1 + 1/cosh[7c(d ffd + X SAT )/(4H sub )]j / 
[1 - 1/cosh [*(d ffd - X SAT )/(4H sub )]] } 



= l ds /W g 



[0103] When the device is dissipating a lot of power and thus unable to outflow 
this heat effectively because of finite thermal conductance, then the device enters a 
thermal regime known as "self-heating". In self-heating, areas adjacent to the device's 
heat sources, or that region approximately under the gate, are heated to higher 
temperatures than the base-plate. As a result, the device becomes hotter than its 
surrounding environment. In this regime, "channel temperature" must be used to gauge 
how hot the intrinsic device gets. 

[0104] As an example of how the formulation accurately models self-heating 

effects, a large periphery 8-fingered 600^m total gate periphery device cell was tested 
for DC-IV and S-parameter at room temperature. This particular device layout 
dissipates too much DC power at high drain voltage or high drain current that could be 
effectively removed by thermal conduction. As a result, the device suffers from 
dramatic self-heating. As shown in Fig. 4, the electrothermal model in the semiphysical 
device model is able to relatively closely simulate the effect of the self-heating channel 
temperature upon the electrical performance of the device. 
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SEMI PHYSICAL MODEL 
[0105] Semi-physical device modeling represents both the physical device 

characteristics and measured characteristics, which can be used to simulate RF 
performance through physically-based device models. The semi-physical model is an 
analytical model based upon empirical expressions that model the physics of HEMT 
operation, hence the terminology "semi-physical". The model incorporates real process 
parameters, such as gate length recess etch depth, recess undercut dimensions, 
passivation nitrite thickness, and the like. By using empirical expressions, the semi- 
physical model is able to maintain relatively good measured to model accuracy while 
accounting for the effects of process variations on the device performance. 
[01 06] The semi-physical model provides model elements for the standard small 
signal equivalent circuit model or FET is illustrated in FIG. 5. FIG. 6 is a rough 
translation of the physical origins for each of the equivalent circuit elements in the 
small-signal equivalent circuit model illustrated in FIG. 5. FIG. 7 is a cross-sectional 
drawing of an exemplary HEMT device structure. However, unlike conventional 
methods, the model elements are derived from small signal excitation analysis of the 
intrinsic charge and electric fields within the device. As such, the simulated small 
signal model elements represent a relatively accurate physical equivalent circuit 
description of a physical FET. 

[0107] The general methodology for the semi-physical modeling of intrinsic 
charge, electrical conductance and the electrical field are as set forth below. First, the 
relationships between the conduction band offsets, electrical permitivities and material 
composition for the various materials in the epi stack are determined. These 
relationships can be performed analytically or by fitting simulated data from physical 
simulators. Subsequently, the basic electron transport characteristics in any of the 
applicable bulk materials in the epi stack are determined. Once the electron transport 
characteristics are determined, the undeleted linear channel mobility is determined 
either through material characterization or physical simulation. Subsequently, the 
Schottky barrier height value or expressions are determined. Once the Schottky barrier 
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height value is determined, the semi-physical equations are constructed modeling the 
following characteristics: 

[01 08] Fundamental-charge control physics for sheet charge in the active channel 
as controlled by the gate terminal voltage. 

[0109] Average centroid position of the sheet charge within the active channel 
width. 

[01 1 0] Position of charge partitioning boundaries as a function of gate, drain and 
source terminal voltages. 

[0111] Bias dependence of linear channel mobility and surface depleted region. 
[0112] Bias dependence of the velocity saturating electric field in the channel 
[0113] Saturated electron velocity. 

[0114] Electrical conductance within the linear region of the channel under the 
gate. 

[0115] Electrical conductance within the source and drain access regions 
[01 1 6] Once the semi-physical equations are determined, the empirical terms of 
the semi-physical modeling equations are adjusted to fit the model I- V (current- voltage) 
characteristics against measured values. Subsequently, the empirical terms are 
interactively readjusted to achieve a simultaneous fit of measured C-V (capacitance- 
voltage) and I-V characteristics. Lastly, the empirical modeling terms are fixed for 
future use. 

[0117] By constructing a comprehensive set of semi-physical equations that cover 

all of the physical phenomenon as mentioned above, the physical operating mechanisms 
within a HEMT device can be relatively accurately determined. FIG. 8 illustrates a set 
of relatively accurate measured-to-modeled I-V characteristics for a HEMT using the 
semi-physical modeling discussed herein. In particular, FIG. 8 illustrates the drain-to- 
source current 1^ as a function of the drain-to-source voltage V ds for various gate biases, 
for example, from 0.4V to -1 .0V. As shown in FIG. 8, solid lines are used to represent 
the semi-physical model while the Xs are used to represent measured values. As shown 
in FIG. 8, a close relationship exists between the measured values and the modeled 
parameters. 
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[01 18] An example of semi-physical modeling for physical device operation in 

accordance with the present invention is provided below. The example utilizes an 
exemplary device as illustrated in FIGS. 9 and 10. Table 2 represents exemplary values 
for the physical cross-section parameters in the model. FIG. 1 1 relates to a blown up 
T-gate characteristic which is correlated to the parameters identified in Table 1 . 
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Table 1 

Values for the Physical Dimension Parameters Input into Device Cross Section 



Layout 


Parameter 


Units 


Value 


Gate Length 


Lg 


[urn] 


0.150 


Wing Length 


Lgw 


mm] 


0.520 


Gate Mushroom Crown Length 


Lgmcl 


[Wn] 


0.200 


Total Gate Height 


Hg 


M 


0.650 


Gate Stem Height 


Hgstem 


[TO] 


0.300 


Gate Sag Height 


Hgsag 


t^rn] 


0.100 


Gate Cross-Sectional Area 


GateArea 




0.187 


Max Cross-Sectional Area 


MaxArea 


[Urn 2 ] 


0.364 


Total Gate Periph 


wg 


[M^mj 


duu.uuu 


# Fingers 


N 


[3 


4.000 


Source-Drain Spacing 


Dsd 


mm] 


1.800 


Gate-Source Spacing 


Dsg 


mm] 


0.700 


Gate-Drain Spacing 


Dgd 


mm] 


1.100 


Gate-Source Recess 


RECsg 


[urn] 


0.160 


Gate-Drain Recess 


RECgd 


mm] 


0.240 


Recess Etch Depth 


Hrec 


[A] 


780.000 


SiN Thickness 


Hsin 


[A] 


750.000 


Gatefeed-Mesa Spacing 


Dgfm 


mm] 


2.000 


Gateend-Mesa Overlap 


Dgem 


mm] 


2.000 


Finger-Finger Spacing Thru Drain 


Dffd 


[Urn] 


16.500 


Finger-Finger Spacing Thru Source 


Dffs 


mm] 


13.500 


Source Airbridge Inset? 


AB? 


[3 


P 


Source Airbridge Inset 


Dsabin 


mm] 


28.000 


Source Airbridge Height 


Hsab 


[urn] 


3.500 


Source-Gate Airbridge Clearance 


Hgsab 


mm] 


1.640 


Source Pad Width 


Ws 


mm] 


12.000 


Drain Pad Width 


Wd 


[Urn] 


14.000 


Substrate Thickness 


Hsub 


mm] 


100.000 



[01 19] As mentioned above, the semi-physical modeling of the intrinsic charge 
and electric field within the HEMT device is initiated by determining the relationships 
between the conduction band offset, electric permitivities and material composition for 
the various materials in the epi stack. Material composition related band offset and 
electric permitivity relationships may be obtained from various references, such as 
"Physics of Semiconductor Devices," by Michael Shur, Prentice Hall, Englewood 
Cliffs, New Jersey 1990. The basic electron transport characteristics, for example, for 
the linear mobility of electron carriers in the bulk GaAs cap layer may be determined 
to be 1350em 2 /Vs, available from "Physics of Semiconductor Devices 55 , supra. The 
linear mobility of electron carriers in the undeleted channels is assumed to be 
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5500cm 2 /V s. This value maybe measured by Hall effect samples whichhave epi stacks 
grown identically to the stack in the example, except for some differences in the GaAs 
cap layer. The Schottky barrier height is assumed to be 1 .05 1 volts, which is typical of 
platinum metal on a AlGaAs material. 

[0120] The following equations represent the semi-physical analytical 
expressions to model the charge control and centroid position in the sample. 



Empirical Charge Control Expression 




[cm' 2 ] 


[1+(N s VN max )v)} 1 ^ 


Ideal Charge Control with Filling Law 


n; 


[cm' 2 ] 


= 2N 0 ln [1+exp(V gt /0iV ttl ))] 


Ideal Charge Control 


No 


[cm 2 ] 


=8! ^ V th / [2 q (d, + A d,) 10000] 






[cm* 2 ] 


=(N ra „ 0 + N 1MxL V (te nnm » 


Maximum Channel Charge 


^max 


Hchan / H chanREF ) 


Initial Gate-Channel Voltage 


V g. 


[V] 


=V ss -O b -AE c -V TO -aV ds 


Threshold Voltage 


V TO 


m 


= %-&E c -V T 


Doping Threshold Voltage 


v T 


m 


= qN sd6lta d 5 10000 1 8, 

- {( H space + H bar + ^fdope + M cap )-H rac 


Gate-to-Channel Spacing 




Em] 


(10 10 ) 




note that the expression for di can be changed for different epi-stacks 








- H ehan [ 1 - d|K*V g te/H chanREF - 


Movement of Sheet Carrier Centroid 




[m] 


^iL*Vd S /H chanREF ] 


Empirical Charge Control Shaping Parameter 


Y 


n 




Semi-Physical Subthreshold Populating Rate 




[] 




Dielectric Permittvity of the Barrier Layer 


Ef 






The thermal voltage 


Vti, 


m 


- K B T amb / q 


Ambient Temperature 


"^"amb 


EK] 




Fixed Emprical Maximum Sheet Charge 


N maxO 


[cm- 2 ] 




Vds Dependent Emprical Maximum Sheet Charge 


^maxL 


[cm" 2 ] 




Vds Dependent Emprical Nmax shaping term 


n Nmax 


II 




Channel Layer Thickness 


^chan 


[A] 




Reference Channel Layer Thickness 


^chanREF 


[A] 






(Channel Thickness for the sample for which the model was first derived) 


Schottky Barrier Height 


Ob 


[V] 




Conduction Band Offset between Channel and Barrier 


A E c 


[V] 




Front Delta Doping 


Nsdelta 


[cm" 2 ] 






note that this expression can be modified for non-delta doped epi -stacks 


Gate-to-Front Delta Doping Spacing 


d8 


[m] 


= {(H bar +H fdope+ H cap )-H rec }/(l0 1c 


Barrier Thickness between front doping and channel 


"space 


[A] 




Barrier Layer Thickness before front doping layer 




[A] 




Front Doping layer thickness 


Hfdope 


[A] 




Cap layer thickness 


H cap 


[A] 




Empirical Drain-Induced Barrier-Lowering Term 


CT 


E3 




Sheet Charge Position Gate Bias Factor 


d iK 


[A/V] 




Sheet Charge Position Drain Bias Factor 


d iL 


[AAV] 








= V th [1+V gl /2V th + 


Effective Gate Voltage 


V gte 


m 


sqrt(o 2 + {V flt /2V th -1) 2 ] 


Empirical Transition Width Parameter 


5 


[] 
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[0121] As used herein, Ns represents the model sheet carrier concentration within 
the active channel. Ns 5 represents the ideal charge control law and is modeled as a 
semi-physical representative of the actual density of state filling rate for energy states 
within the channel v. gate voltage. The gate-to-channel voltage used for the charge 
control, Vgt, is a function of the Schottky barrier height, conduction band offsets and 
doping in the epi stack as is known in the art. 

[0122] The following equations represent the semi-physical expression used to 
model the position of regional charge boundaries within the HEMT device. These 
expressions govern how to partition the model charge between the influence of different 
terminals. 



w 


Effective Gate Length 


L geff 




= L g + A L B + A L d 




Gate-Source Control Region 






= L B /2+ A L s + 




Source-Side Effective Gate Length Extension 


Al-. 




- aL s o + aLk*V 




Dram-Side Effective Gate Length Extension 




[(i»"J 


= A«- d0 + aL k * V gte + A L L * V dse2 




Gate-Drain Control Region 


L 0d 


t(i»n] 


= (L g /2+ A L d )*{tanh[10{L g /2-X D1 )] + 1}/2 










= X^V^M^/ 










{M„ K V ata *(1+IX 0L V lh M )ML / 




Bias Dependent Extension of the Saturated Transport Region 


^D1 


EjimJ 






Empincal Drain-Saturated Transport Boundary Factor 






= L g V ds / { 2 [ 1 + ]< 1/m > } 




















V ds m (V ds /V 3Htn )^ 1 V{ 










V salfl m[1 + (V d A„ tr ,n (1/ra+,) }] 




Position of the Boundary between Regions 1 and 2 




Efini] 


V ds Mx sL + V gte M X5K } 




Note: Region 1 denotes the linear region, while Region 2 c 


le notes the saturated region of the channel 




Empirical Effective Gate Length Extension Gate Bias Factor 










Empincal Effective Gate Length Extension Drain Bias Factor 


Ak 


E^m/V] 






Effective Drain-Source Voltage Control-2 


Vdse2 


M 


-v ds /[i+(v d5 /v satn n i/m 




Rough, Intrinsic Saturation Voltage 


Vsatn 


rvj 


- ! sat^9ehl 
= 9ehl Vgte / 




Rough, Intnnsic Saturation Current Level 




[A] 


[1 + g^R, + sqrt(1+2g clu ^+(V B(B A/ L ) 2 )I 




Intnnsic Conductance of the Linear Region, Under the gate 


9chi 


[S] 


=(q N s H ave Wg)/Lg 




Rough Intnnsic Saturation Voltage Level 


v L 


[V] 


= F s *L g 




Empincal Knee Shaping Parameter 


m 


El 






Ernpincat Region 2 extension Drain Bias Factor 


M XdL 


n 






Empincai Region 2 extension Gate Bias Factor 


MxdK 


n 






Fine Intnnsic Saturation Voltage 


Vsaten 


IV] 












= 9ch 1 VV L M-V L {A + g ohl R E )H- 










sqrt(V L 2 (A+g 6hl R s ) 2 +V gt9 2 -{g chl R s V L ) 2 )] / 




Fine Intnnsic Saturation Current Level 


'satcom 


[A] 


[V 2 ( 1-g ch ,R s (V L /V gte } 2 )] 




Saturation Region Length Ratio 


A 


u 


= X s /L geff 




Initial Starting position for Region 1 & 2 Boundary 


x s0 




= Lg/2 




Region 1 & 2 Boundary Bias Factor 




n 






Region 1 & 2 Boundary Drain Bias Factor 




EI 






Region 1 & 2 Boundary Drain Bias Factor 


M XsK 


El 
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[0123] The following equations represent is the semi-physical expression used 
to model the bias dependence of linear channel mobility in depleted regions. 



Depleted Channel Mobility u. ave [cm 2 /V*s] = M dchan + I^K-Vgte 

Fixed Depleted Channel Mobility u.dchan l cm2/v * s l 
Depleted Channel Mobility Gate Bias Factor h<jk [cm 2 /V 2 *s] 



[0124] The following equations are the semi-physical expressions used to model 

the bias dependence of the saturating electric field and saturation velocity. 



= v s I 

Saturating Electric Field F s [V/^m] [(n sat + (WW 1 0000] 
Fixed Saturating Channel Mobility u, sat [cm 2 /V*s] 
Saturating Channel Mobility Gate Bias Factor fi sat K [cm 2 A/ 2 *s] 
Saturation Velocity v s [cm/s] 



[0125] FIG. 12 is a schematically illustrates how electrical conductance in the 
source and drain access regions are modeled in the example. 

[0126] The following equations describe the semi-physical model for the source 

access region conductance: 
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Source Access Resistance K s En] 
Source Access Resistance 1 Channel and Cap R s undepCap En*ia m ] 

Source Access Resistance: Recess and Undepleted Cap R SAccess Ifi*p. m J 



Source Access Resistance: Crowding resistance due to conductance mismatch 
Resistance of the Source Recess Access region at high on-state bias (V^) 
Resistance of the Undepleted Source Recess Access region 
Uncapped, Fully Depleted Sheet Resistance 
Capped, Undepleted Sheet Resistance 
Uncapped, Undepleted Sheet Resistance 
Cap Sheet Resistance 
Surface Depletion Factor 
High On-state bias: Diode Turn-on voitage 
Ohmic Contact Resistance 
RF Ohmic Contact Resistance Reduction Factor 
Source Access Resistance Bias Modification Factor 
Cf-Vds Bias Modification Factor 
Rs-Vds Bias Modification Factor 
Rs-Vgs Bias Modification Factor 
Cf-Vgs Swith point to On-state 
Cf-Vgs Bias Expansion Factor 
Rs-Vgs Swith point to On-state 
Rs-Vgs Swith point to Off-state 
Rs-Vds Swith point from Off-On transition 
Rs-Vds Bias Expansion Factor 
Rs-Vgs Bias Expansion Factor 
Rs-Vds Bias Expansion Factor @ Rs Saturation 
Rs Bias Shaping Factor 



" SBaundary 
R ON 
" S depRec 

R S undepRec 

^SHdep 

R SH 

R SHundep 

R SHCap 

F 

1 surfUndep 

R cont 
R^rconF 

MR S 

MR sL 
MR sK 

vc^ 

VR S0 „ 
VR SKnee 

KR sL 
KR sK 

KR sSat 
yR s 



[fi/sq] 
fo/sq] 
[Q/sq] 

Ea-sq] 

n 
m 



~ ( Rsund«pCap + ^SAccess + RsBoundary ) ' ^V g 
= Rcont^RFrconp + 

R SH [ D SEJ -(REC sg +L g /2)] 

= R 8(toP R.e 0N *M^*tanh{tKC«* ( 
- VC f0n ^V d /MC fL )] + 1 } / 2 * { V B3 

1 2 *[1-tanh< KR, K <V g , - VR, 0n )))r { tanh { 

KR, Sal (Vd» " VR SK „«)] + 1 } / 2 + RsundepRec * { fcnh[ 

KR,k(V s ,-VR, 0 „}] + 1}/2 
= Rsd eP Rec° N *MR t *tanh{[KR, K *V BS + KR sL *V d ^ 
VR s0 fr ] + 1 } / 2 * {{1 +V d ,MR, L )*MR, K *[1- 

tanh{KR, Sat (V ds .VR sKn8e )}]/ [ 2*(1 + [V fls /((1 + 

V ds MR fL rMR sK )]^/ 1/ v Rs >]} 

= ^SHeiep ( R^C sg ) 
= RsHundep ( REC sg ) 

- 1 I {q N max ^ ave ) 

= 1/{1/RsHCap + 1/R S Hunde P ) 
= ^surtUndep ^ (M ^max Kjndchan) 

= 1 / [q N sCap H cap ( H cap - H capEleh ) ] 



= <t> b -A Ec -A Ef 



[0127] The following equations describe the drain access. 
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Drain Access Resistance 

Drain Access Resistance, except for recess access 
Drain Recess Access Resistance 



Enl 

Cq] 
Eol 



Drain Access Resistance: Recess and Un depleted Cap R DAccess lfi*u. m l 

Drain Access Resistance 1 Channel and Cap R p undepCap ln*jx m ] 
Resistance of the Undepleted Drain Recess Access region R D undepRec fo*^ m l 



Resistance of the Saturated Drain Recess Access region R ossturated fo*n m 3 

Resistance of the Drain Recess Access region at high on-state bias (V on ) R D de pR ec 0N fo*u™l 

R DRec 

Drain Access Resistance Bias Modification Factor MF^ [ ] 

Rd-Vds Access Bias Modification Factor MR daccL [ ] 

Rd-Vgs Swith point to On-state VR d0n [ ] 

Rd-Vgs Swith point at saturation VR dSat E ] 

Rdaccess-Vgs Bias Expansion Factor KR daccK [ ] 

Rd-Vds Saturation Bias Modification Factor MR dSat E ] 

Rd-Vds Access Resistance Bias Modification Factor MR dL [ } 

Rd-Vgs Saturation Swith point to Off-state VR dSat0ff [ ] 

Rd-Vgs Saturation Swith point VR dSat [ ] 

Rd-Vgs Saturation Swith point transition width VR dSata [ ] 

Rd Bias Shaping Factor yR d E 3 



- ( ^DRF + RdN ) 

! ( Rd Access + RDundepCap + RoundepRec + RDSaturated ) ^ ^V g 

- Rdn 

= ( RDRec) ^ g 
- R D depRec ON (MR d *V 3S -M)[1- 

tanh{ KR daecK * [ V gs - VR dDn + V ds *MR daccL }}]/2 
= R C on/RF rconF + R SH [ D gd 

- (REC gd +L a /2)] 

~ RsHundep ( REC gd ) 

= MBWV^expf -JV^+Vj/MR^+VR^s^] 2 / (2 
} / { VR^^ [1 +(V ds A/Rd Sat )' Rd ]< ^ RdA + 1 '* Rd) } 

= RsHdep { REO gd ) 
= SVL ( 1 + ( RDAcce„ ' RDundepCap V( 2Svl) + S* 

+ ((RdAccsss " RDundepCap y(26vi.) ""!)]} 



SEMI-PHYSICAL DETERMINATION OF SMALL-SIGNAL EQUIVALENT 

CIRCUITS 

[0128] To derive values for the familiar small signal equivalent circuit as shown 
in FIG. 5, a small signal excitation analysis must be applied to the semi-physically 
modeled physical expressions. The method of applying such an analysis is as follows: 
[0129] 1) Gate Terminal Voltage Excitation 

[0130] a) Apply a small+Z- voltage delta around the desired bias 

condition, across the gate-source terminals. 

[0131] b) Equivalent circuit element Gm = delta(Ids)/delta (Vgs') where 

delta (Vgs') is mostly the applied voltage deltas, but also subtracting out that voltage 
which is dropped across the gate source access region, shown as RsCont, RsundepCap, 
RsundepRec, ResdepRec, and RsBoundary in FIG. 12, above. 

[0132] c) Equivalent circuit element Cgs and Cgd takes the form of 

delta(Nsn)/delta(Vgs)*Lgn, where delta (Nsn) is the appropriate charge control 



23 



expression, and Lgn is the gate source or gate drain charge partitioning boundary 
length. 

[0133] d) Equivalent circuit element Ri=Lgs/(Cgschannel * vs) where 

Cgs channel is the portion of gate source capacitance attributed to the 
channel only, and vs is the saturated electron velocity. 

[0134] 2) Drain Terminal Voltage Excitation 

[0135] a) Apply a small +/- voltage delta around te same bias condition 

as in 1, but the delta is applied across drain source terminals, 

[0136] b) Equivalent circuit element Rds = l/{delta(Ids)/delta(Vds')} 

where Vds' is mostly the applied voltage deltas, but also subtracting out voltage which 
is dropped over both the gate source and gate drain access regions. 
[0137] c) Equivalent circuit element Cds is taken to be the sum of the 

appropriate fringing capacitance Semi-Physical models, or can take the form of 
delta(Nsd)/delta(Vds 5 )*Xsat, were Nsd is the charge control expression for charge 
accumulation between the appropriate source and drain charge boundaries, and Xsat is 
the length of the saturated region, if in saturation. 

[0138] 3) On-mesa Parasitic Elements: The equivalent circuit elements, Rs and 
Rd are expressed by the appropriate electrical conduction models of the source and 
drain access regions. 

[0139] The RF performance can be predicted at an arbitrary bias point. 
[0140] Table 3 represents a comparison of the values for a high frequency 
equivalent circuit model derived from equivalent circuit model extraction from and 
semi-physical modeling for the sample illustrated in Table 2. 
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Table 2 

Comparison of Modeled Equivalent Circuit Results for Semi-physical 
Modeling Method, and Equivalent Circuit Model Extraction 



Intrinsic Equivalent 
Circuit Parameter 


Equivalent Circuit 
Model 


bemi-rnysicai 
Device Model 


Cqs 


0.227745 pF 


0.182 pF 


Rgs 


64242 Q 


infinite Q 


Cqd 


0.017019 pF 


0.020 pF 


Rgd 


133450 Q 


infinite Q 


Cds 


0.047544 pF 


0.033 pF 


Rds 


160.1791 Q 


178.1 n 


Gm 


135.7568 mS 


124 mS 


Ri 


3.034 Q 


2.553 Q 


Tau 


0.443867 pS 


0.33 pS 



[0141] The results of the semi-physical modeling method produce a small-signal 
equivalent circuit values which are relatively more accurate than the physical device 
simulator in this case. Furthermore, given the differences in the parasitic embedding, 
treatment of the two approaches, the results given in Table 2 yield much closer results 
than a comparison of equivalent circuit values. 

[0142] Table 3 lists the values of parasitic elements used in the model 

derivations. An important difference between the extracted equivalent circuit model 
and the semi-physically derived one is the use of Cpg and Cpd to model the effect of 
launch capacitance for the tested structure. This difference leads to the results of the 
extracted model results being slightly off from the optimum physically significant 
solution. 

Table 3 

Comparison of Modeled "Parasitic" Equivalent Circuit Results for Semi- 
physical Modeling Method, and Equivalent Circuit Model Extraction 



Extrinsic Equivalent 
Circuit Parameter 


Equivalent 
Circuit Model 


Semi-Physical 
Device Model 


Rg 


1 .678 Q 


ua 


Lfl 


0.029314 nH 


0.03 nH 


Rs 


1.7 Q 


1.21 Q 


Ls 


0.002104 nH 


0.003 nH 


Rd 


3.309899 Q 


5.07 n 
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Extrinsic Equivalent 
Circuit Parameter 


Equivalent 
Circuit Model 


Semi-Physical 
Device Model 


Ld 


0.031671 nH 


0.02 nH 


Cpg 


OpF 


0.02 pF 


Cod 


0 pF 


0.01 pF 



[0143] As shown in FIGS. 13, 14 and 15, the modeled results that are simulated 
using the semi-phy sically derived equivalent circuit model very accurately replicate the 
measured high frequency, S-parameter data. 

[01 44] The following equations represent the small-signal excitation derivation 
of small-signal equivalent circuit modeled Gm. Figure 56 illustrates the semi- 
physically simulated bias equations of the small signal Gm compared to measured date. 



=gcnV,, s (mv ds )/ 



Semi-Physicaily Modeled Drain-Source Current Control 




[A] 


[i+(v ds /v sate n i/m 


Small-Signal Deterimination of equiv. Circuit Gm value 


9mRF 


[S] 


= c/l ds /d(V gs -V sAcc ) 






- *dsW* ( R SUndepCap + R SAccess + R S 


Source-Access voltage drop 


^SAcc 


[VI 


RprobeS^g ) 


Fine Extrinsic Saturation Voltage 


V sate 


[V] 


= Uatcom^ 9ch 








QchJ 


Extrinsic Conductance of the Linear Region, Under the gate 


9 ch 


[S] 


[i+gchi(R s + Rd>3 


Intrinsic Conductance of the Linear Region, Under the gate 


9chi 


[S3 





[0145] The following equations represent the small-signal excitation derivation 

of Rds. Figure 57 illustrates the semi-physically simulated bias- 
dependence of the small-signal Rds. 



Small-Signal Deterimination of equiv. Circuit Rds value 


Rd, 


In] 




9dsRF 


[S3 


Drain-Access voltage drop 


y 

v DAcc 


m 


Drain-Saturated Region voltage drop 


^DSAT 


m 


External Test probe or lead resistance 


P 

^probeD 


In! 


High Frequency conductance dispersion factor 


r dsF 


[] 


High Frequency conductance dispersion 


RF rdsF 


[] 



= 1 / 9dsRF 

= { Cf l d5 / d ( V ds - RprobeD * ! ds * V sAcc- V dAcc- V dSat ) 

> * r dsF 

- ^dsW * ( ^DUndepCap + ^DUndepRec * Access + 
Rprobec/^g ) 

^Saturated ) 



= (RFrd*F + 1)* 

tanh(10*iV d5 -V th |) + 1 
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[0146] 



The following equations may be used for illustrating small-signal 
excitation derivation of Cgs and Cgd. Figure 58 illustrates the Semi- 
Physically simulated bias-dependence of the small-signal Cgs and Cgd. 



Small-Signal Deterimination of equiv. Circuit Cgs value 








c <* 


[fF/um] 


Parasitic Gate-Source Fringing Capacitance 


C 8 sf 


[fF^m] 


Parasitic Gate-Source Fringing Capacitance 




[fF/^m] 


Total Specific Gate-Channel Capacitance 


^gcTot 


[fFyn 2 J 


Specific Gate-Channel Capacitance 


C BC 


[fF/^m ] 


Effective Drain-Source Voltage Control 


Vdse 


m 


Specific Gate-Donor Layer Accumulation Capacitance 


^ gcdonor 


[fF/^m 2 ] 


ideal Specific Gate-Channe! Capacitance 


C 9C 


[fFVm 2 ] 


Empirical Parasitic Donor Charge Control Expression 


N 

sDonor 


[cm" 2 ] 


Fringing capacitance to surface of source-access region 


C gswf 


[fF/^m] 


Empirical Fringing capacitance-bias shaping expression 


Cf-form1 


E3 


Fringing capacitance to source-access region 


^ gsfSource 


[fF/^m] 


Fringing capacitance to source metal pads 


^ gsfPad 


[fF/^m] 


Fringing capacitance to drain-access Capped Region 


^ gdfCap 


[fF/^m] 


Fringing capacitance to drain metal pads 


GgdfPad 


[fF/^mj 


Dielectric Coating Thickness Factor 


SiNF 


tl 


Specific Gate-Channel Capacitance Bias ModiflcationFactor 


Msh Kcahn 


[] 


Donor Charge Bias ModiflcationFactor 


M sh Kpar 


[] 


Empirical Specific Charge Control Shaping Parameter 


Y c 


[] 



{ 1 " [ (V sat en-V dse ) / ( 2*V satan -V dse )f } 



= c C( 



{1-[V satan /{2*V satBn -V dse )] 2 } 
= Cg surf C f . form1 SiNF + Cgsf Source + Cgsf Pad 
= Og surf C f . form1 SiNF + Cgdf Cap + Cgdf Pad 

— Cg C + Cg cc | onor 

c gc Msh Kchan / 

= V ds /[1 + (V ds /V saten n 1/ - 
= qdN sDonor /dV gs 
= q dN s '/cfV gs 
KN.' + N.m + AdjMsh^^V^/ 
(d 1 + A dl ) 



{ 1 - tanh[ KC, K ( V gs - VC fQn + V ds MC fL )] } / 2 



[0147] The following equations are involved in the small-signal excitation 
derivation of Ri Figure 59 which follows shows the Semi-Physically 
simulated bias-dependence of the small-signal Ri. 



Gate-Source Non-quasistatic charging resistance ^icharge En*ja m 3 - L gs2 W g / [ C gsChan W g v s ] 

~ C gsf p C gc L gs * 

Gate-Channel Source Capacitance C gsChan [fF^m] { 1 - [ (V saten -V dse ) / ( 2*V saten -V dse )] 2 } 
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EXAMPLE OF SEMI-PHYSICAL MODEL AND BIAS DEPENDENCE 



SMALL SIGNAL SOURCE AND DRAIN RESISTANCE. RS AND RD 
[0148] FIG. 60 shows the semi-physically simulated bias dependence of the on- 

mesa parasitic access resistances Rs and Rd. 

[01 49] The following example verifies how the semi-physical small-signal device 
model is able to provide accurate projections for bias-dependent small-signal 
performance. In this example, the same semi-physical device model as used in the 
previous examples was used because the example MMIC circuit was fabricated utilizing 
the same HEMT device technology. 

[0150] In this example, the bias-dependence small-signal gain and noise 
performance of a two-stage balanced K-band MMIC LNA amplifier is replicated 
through microwave circuit simulation using small signal and noise equivalent circuits 
that were generated by the semi-physical model. The results of the measured and 
modeled results are shown below in Table 4. As seen from these results, the Semi- 
Physical device model was able to accurately simulate the measured bias-dependent 
performance, even though the bias variation was quite wide. 



Table 4 

Measured vs Modeled NF and Gain @ 23.5 GHz for a K-band MMIC LNA at 

different bias conditions. 



Bias Condition 


Measured Gain 


Predicted Gain 


Measured NF 


Predicted 




@ 23.5 GHz 


@ 23.5 GHz 


@ 23.5 GHz 


NF @ 23.5 
GHz 


Vds=0.5 V 


15.2 dB 


15.8 dB 


2.97 dB 


2.77 dB 


112 mA/mm 










Vds=1.0 V 


20.6 dB 


21.0 dB 


2.29 dB 


2.20 dB 


112 mA/mm 










Vds=2.0 V 


19.8 dB 


20.2 dB 


2.25 dB 


2.15 dB 


112 mA/mm 










Vds=3.0 V 


18.9 dB 


19.1 dB 


2.30 dB 


2.11 dB 


112 mA/mm 










Vds=3.5 V 


18.4 dB 


18.5 dB 


2.34 dB 


2.18 dB 
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112 mA/mm 










Vds=4.0 V 


18.0 dB 


18.0 dB 


2.37 dB 


2.27 dB 


112 mA/mm 










Vds=2.0 V 


16.4 dB 


18.0 dB 


2.45 dB 


2.21 dB 


56 mA/mm 










Vds=2.0 V 


21.4 dB 


20.9 dB 


2.38 dB 


2.21 dB 


170 mA/mm 










Vds=2.0 V 


22.2 dB 


21.0 dB 


2.65 dB 


2.6 dB 


225 mA/mm 










Vds=3.0 V 


21.4 dB 


20.3 dB 


2.71 dB 


2.61 dB 


225 mA/mm 










Vds=3.0 V 


20.5 dB 


20.0 dB 


2.42 dB 


2.22 dB 


170 mA/mm 










Vds=4.0 V 


19.6 dB 


19.2 dB 


2.50 dB 


2.29 dB 


170 mA/mm 











[0151] A plot of measured vs modeled gain for the values listed in Table 3 above 
is shown in Figure 61 . 

[0152] The following example verifies how the Semi-Physical small-signal 
device model is able to provide accurate projections for physically dependent small- 
signal performance. In this example, the same Semi-Physical device model as used in 
the previous examples was used. 

[0153] In this example, physical process variation was input into the Semi- 
Physical device model in terms of statistical variation about known averages, cross- 
correlation, and standard deviations. The goal of this exercise was to replicate the 
measured DC and small-signal device variation. The degree of accurate replication 
indicates the degree to which the Semi-Physical model is physically accurate. 
[0154] Table 5 below lists the simulated, and known process variation that was 
used: 
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TABLES 
Statistical process variation model 

Parameter Nominal Standard Dev. 

Gate Length 0,15 um 0.01 urn 

Gate-Source Recess 0.16 um 0.015 um 

Gate-Drain Recess 0.24 um 0.020 um 

Etch Depth 780 A 25 A 

Pass. Nitride Thickness 750 A 25 A 

Gate-Source Spacing 0.7 um 0.1 um 

Source-Drain Spacing 1.8 um 0.15 um 



[0155] In the course of microelectronic component production, sample devices 
are tested in process in order to gain statistical process control monitor (PCM) data. 
Figure 1 8 shows schematically the kind of data that is extracted and recorded from 
measured device I-V's during PCM testing. 

[0156] Since the Semi-Physical device model is able to simulate I-V's, it was 
able to simulate to variation of I-V's due to physical process variation. These I-V's 
were analyzed in the same fashion to extract the same parameters that are recorded for 
PCM testing. Figures 63, 64 AND 65 shows how accurately the simulated results 
match with measured process variation. Figure 63 shows how the Semi-Physically 
simulated Vgpk and Gmpk match with actual production measurements. Figure 64 
shows how simulated Idpk and Gmpk match, also. Finally, Figure 65 shows how 
simulated Imax and Vpo also match very well. 

[0157] Small-signal S-parameter measurements are also taken in process for 
process control monitoring. These measurements are used to extract simple equivalent 
circuit models that fit the measured S-parameters. Since the Semi-Physical device 
model is able to simulate these equivalent circuit models, it was able to simulate the 
variation of model parameters due to physical process variation. 
[0158] Figures 66 and 67 show how accurately the simulated results match with 
measured/extracted process variation for the small-signal model parameters. Figure 
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66 shows how the Semi-Physically simulated Rds and Gm match very well with actual 
extracted model process variation. 

[0159] More direct and convincing evidence supporting the accurate, physical 
nature of the Semi-Physical model can be shown be comparing the dependence of 
simulated and measured performance to real physical variable. As shown in Figure 68 , 
the Semi-Physical model is able to very accurately reproduce the dependence of Imax 
upon gate length. In addition, the Semi-Physical model is also able to replicate 
physical dependence for high-frequency small-signal equivalent circuits. This is 
shown in Figure 69, which shows that it is able to reproduce the dependence of Rds 
with Recess undercut width. 

S-PARAMETER MICROSCOPY 
[0160] The S-parameter (SPM) method utilizes bias dependent S-parameter 
measurements as a form of microscopy to provide qualitative analysis of the internal 
charge and electrical field structure of the semiconductor device heretofore unknown. 
Pseudo images are gathered in the form of S-parameter measurements extracted as 
small signal models to form charge control maps. Although finite element device 
simulations have heretofore been used to calculate the internal charge/electric field of 
semiconductor devices, such methods are known to be relatively inaccurate. The S- 
parameter microscopy provides a relatively accurate method for determining the 
internal charge and electric field within a semiconductor device. With accurate 
modeling of the internal charge and electric field, all of the external electrical 
characteristics of semiconductor devices can be relatively accurately modeled 
including its high frequency performance. Thus, the system is suitable for making 
device technology models that enabled high frequency MMIC yield analysis 
forecasting and design for manufacturing analysis. 

[0161] S-parameter microscopy is similar to other microscopy techniques in that 
SPM utilizes measurements of energy reflected to and from a sample to derive 
information. More particularly, SPM is based on transmitted and reflective microwave 
and millimeter wave electromagnetic power or S -parameters. As such, S-parameter 
microscopy is analogous to the combined operation of scanning and transmission 
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electron microscopes (SEM and TEM). Scattered RP energy is analogous to the 
reflection and transmission of the electron beams in the SEM and TEMs. However, 
instead of using electron detectors as in the SEM and TEMs, reflectometers in a 
network analyzer are used in S-parameter microscopy to measure a signal. S- 
parameter microscopy is similar to other microscopy techniques in that both utilize; 
measurement of scattering phenomenon as data; include mechanisms to focus 
measurements for better resolution; and include mechanisms to contrast portions of the 
measurement to discriminate detail as shown in Table 6 below: 



Table 6 



General Microscopes 


S-Parameter Microscope 


Measure of scattered energy 


Measures S-Parameters 


Mechanism for "focus " 


Focuses bv extraction of Unique 
equivalent circuit models 


Mechanism for "contrast" 


Contrasts by using bias dependence to 
finely discriminate the nature and 
location of charge/electric fields 



[0162] RESULT : Detailed "images" of device's internal charge and electric 
field structure. 

[0163] Images as discussed herein, in connection with S-parameter microscopy 
do not relate to real images, but are used provide insight and qualitative detail 
regarding the internal operation of a device. More specifically, S-parameter 
microscopy does not provide visual images as in the case of traditional forms of 
microscopy. Rather, S-parameter microscopy images are more like maps which are 
computed and based on a non-intuitive set of measurements. 

[0164] FIG. 16 illustrates a conceptual representation of an S-parameter 
microscope, generally identified with the reference numeral 20. The S-parameter 
microscope 20 is analogous to a microscope which combines the principles of SEM 
and TEM. Whereas SEM measures reflections and TEM measures transmissions, the 
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2-port S-parameter microscope 20 measures both reflective and transmitted power. As 
a result, data derived from the 2-port S-parameter microscope contains information 
about the intrinsic and extrinsic charge structure of a device. More particularly, as in 
known in the art, SEM provides relatively detailed images of the surface of a sample 
through reflected electrons while TEM provides images of the internal structure 
through transmitted electrons. The reflective signals are used to form the external 
details of a sample while transmitted electrons provide information about the interior 
structure of a device. In accordance with an important aspect of the invention, S- 
parameter microscopy utilizes a process of measuring reflective and transmitted signals 
to provide similar "images" of the charge structure of a semiconductor device. As used 
herein the internal and external electrical structure of a semiconductor device are 
commonly referred to as intrinsic device region 22 and extrinsic parasitic access region 
24 as shown in FIG. 17. Also contributing to the external electrical structure of the 
device are parasitic components associated with the electrodes and interconnects which 
are not shown. These are the so-called device "layout parasitics" 
[0165] Referring to FIG. 16, the ports 26 and 28 are emulated by S-parameter 
measurements. The S-parameter measurements for a specific semiconductor device, 
generally identified with the reference number 30, are processed to provide charge 
control maps, shown within the circle 32, analogous to images in other microscopy 
techniques. These charge control maps 32, as will be discussed in more detail below, 
are expressed in the form of equivalent circuit models. As shown in FIG. 18, linear 
circuit elements are used in the models to represent the magnitude and state of 
charge/electric fields inside the semiconductor device 30 or its so-called internal 
electrical structure. The position of the circuit elements within the model topology is 
roughly approximate the physical location within the device structure, hence the charge 
control map represents a diagram of the device's internal electrical structure. 
[0166] The interpretation of the exact location of measured charge/electric fields 
within the semiconductor device is known to be ambiguous since an equivalent circuit 
model, for example, as illustrated in FIG. 19 with discrete linear elements, is used to 
represent the distributed structure of the charge/electric fields in the actual device. 
Although there is no exact method for distinguishing the physical boundaries between 
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measured quantities, bias dependence is used to clarify how the S-parameters should 
be discriminated, separated and contrasted. In particular, changing bias conditions is 
known to change the magnitude and shift boundaries between the charge and electric 
fields within the device. The changes are normally predictable and qualitatively well 
known in most technologies. As such, the charge control maps can readily be used as 
maps illustrating the characterization of physical changes in magnitude, location and 
separation of electric charge and electric fields. 

[01 67] Analogous to other forms of microscopy, the S-parameter microscope 20 
in accordance with the present invention also emulates a lens, identified with the 
reference numeral 40 (FIG. 16). The lens 40 is simulated by a method for the 
extraction of a unique equivalent circuit model that also accurately simulates the 
measured S-parameter. More particularly, parameter extraction methods for equivalent 
circuit models that simulate S-parameters are relatively well known. However, when 
the only goal is accurately fitting measuring S-parameters, an infinite number of 
solutions exist for possible equivalent circuit parameter values. Thus, in accordance 
with an important aspect of the present invention, only a single unique solution is 
extracted which accurately describes the physical charge control map of the device. 
This method for unique extraction of equivalent circuit model parameters acts as a lens 
for focus the charge control map solution. As discussed and illustrated herein, the lens 
40 is subsequently simulated by a filter that is based on an apparent layout parasitic 
embedding model. As discussed below, the layout parasitic embedding model consists 
of linear elements which simulate the effect of the device's electrodes and 
interconnects upon its external electrical characteristics. A Pi FET embedding model 
42, is described below. This model effectively acts as a filter to remove the electrical 
structure of the extrinsic parasitic access contribution to the preliminary charge control 
map solution. The resultant, filtered charge control map solution represents a clearer 
"image", which shows only the electrical structure of the intrinsic device. This 
enhanced imaging is needed in order to achieve as accurate a view of the internal 
electric charge/field as possible. Unlike conventional extraction techniques as 
illustrated in FIG. 2 1 , which can only extract equivalent non-unique circuit models and 
not the unique charge control map, the S-parameter microscope 20 in accordance with 
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the present invention is able to relatively accurately model the internal electric 
charge/field structure within a semiconductor device. 

[01 68] An exemplary application of the S-parameter microscope is illustrated in 
detail below. In this example, an exemplary GaAs HEMT device with four gate 
fingers and 200jim total gate periphery formed in a Pi-FET layout as generally 
illustrated in FIG. 22 and identified with the reference numeral 43 is used. The GaAs 
HEMT 43 is adapted to be embedded in a 100-|im pitch coplanar test structure to 
facilitate on wafer S-parameter measurement. 

[0169] Initially, as illustrated in FIGS. 23 and 24, the I-V characteristics for the 
device are measured. In particular, the drain source current Ids is plotted as a function 
of drain-to-source voltage Vds at various gate voltages Vgs as shown in FIG. 23. and 
FIG. 24 illustrates the drain-to-source current Ids as a function of gate voltage Vgs and 
transconductance Gm (i.e. the derivative of Ids with respect to Vgs) at different drain 
voltages Vds. These I-V characteristics are typical of HEMT devices and most 
semiconductor devices, which are one type of three-terminal semiconductor device 
technology. 

[01 70] Table 7 shows the bias conditions in which S-parameters were measured. 
The S-parameters were measured from 0.05 to 40 GHz at each bias condition. FIG. 
25 illustrates a Smith chart illustrating the measured S-parameters SI 1, S 12 and S22 
for frequencies from 0.05 to 40.0 GHz. FIG. 26 is a graphical illustration of magnitude 
as a function of angles for the measured S-parameter S21 for frequencies from 40.05 
to 40.0 GHz. 
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TABLE 7 

Measured S-parameter Bias Conditions 



Biases 


Vds = 
0 V 


Vds = 
0.5 V 


Vds = 
1.0 V 


Vds = 
2.0 V 


Yds = 
4.0 V 


Vds = 
5.0 V 


Vgs 


-1.6 V 


Yes 


Yes 


Yes 


Yes 


Yes 


Yes 


-1.4 V 


Yes 


Yes 


Yes 


Yes 


Yes 


Yes 


-1.2 V 


Yes 


Yes 


Yes 


Yes 


Yes 


Yes 


-1 V 


Yes 


Yes 


Yes 


Yes 


Yes 


Yes 


-0.8 V 


Yes 


Yes 


Yes 


Yes 


Yes 


Yes 


-0.6 V 


Yes 


Yes 


Yes 


Yes 


Yes 


Yes 


-0.4 V 


Yes 


Yes 


Yes 


Yes 


Yes 


Yes 


-0.2 V 


Yes 


Yes 


Yes 


Yes 


Yes 


Yes 


0 V 


Yes 


Yes 


Yes 


Yes 


Yes 


Yes 


0.2 V 


Yes 


Yes 


Yes 


Yes 


Yes 


Yes 


0.4 V 


Yes 


Yes 


Yes 


Yes 


Yes 


Yes 


0.6 V 


Yes 


Yes 


Yes 


Yes 


Yes 


Yes 



[0171] Using the small signal model illustrated in FIG. 19, the extracted small 
signal equivalent circuit values are obtained as illustrated in Table 8 for each S- 
parameter at each bias condition, using the extraction method discussed below. 
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[0172] The values in Table 8 represent solutions that are close to the charge 
control map and represent physically significant solutions of the FET's electrical 
structure. However, the values represented in Table 8 contain the influence of external 
layout parasitics which are subtracted using a model for the embedding parasitics to 
obtain the most accurate charge control mapping to the intrinsic device characteristic. 
In particular, an embedding model is applied filter the extracted equivalent circuit 
model values and obtain values more representative of the intrinsic device. In 
particular, in the exemplary embodiment, aPiFET embedding parasitic model is used 
to subtract capacitive contributions due to interelectrode and off-mesa layout parasitic 
influences. This filter essentially subtracts known quantities formed from the 
parameters Cgs, Cgd and Cds depending on the device layout involved. In this 
example, embedding of the inductive parameters is not necessary because these 
quantities are extrinsic and do not contribute to the charge control map of the intrinsic 
device, 

[01 73] As discussed above, the lens with filter are used to generate unique charge 
control maps. In particular, FIGS. 27-30 illustrate the bias dependent charge control 
maps for the parameters RS, RD, RI, CGS and CGD as a function of bias. More 
particularly, FIG. 27 illustrates a charge control map of the charge and electric field 
distribution in the on-mesa source access region illustrated by the source resistance It, 
as a function of bias. FIG. 28 illustrates a charge control map of the charge and 
electric field distribution in the on-mesa drain access region illustrated by the drain 
resistance R d as a function of bias. FIG. 29 illustrates a charge control map for a non- 
quasistatic majority carrier transport illustrated by the intrinsic device charging 
resistance R { as a function of gate bias for different drain bias points. FIG. 30 
illustrates a charge control map for gate modulated charge and distribution under the 
gate shown with the gate capacitance CGS and CGD as a function of bias. 
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FILTER 

[0174] As mentioned above, the S-parameter microscope 20 may utilize a filter 
to provide a clearer charge control map for modeling the internal electric charge/field 
of a semiconductor device. Although the filter is illustrated in connection with the 
PiFET with multiple gate fingers, as illustrated in FIGS. 31 and 32, the principles of 
the invention are applicable to other semiconductor devices. 

[01 75] As illustrated in FIG. 3 1 , PiFETs are devices in which the gate fingers and 
the edge of the active region resemble the greek letter it, as illustrated. Such PiFET 
layouts facilitate construction of multi fingered large periphery device cells, for 
example, as illustrated in FIG 32. In accordance with an important aspect of the 
invention, the multi-finger semiconductor device is modeled as a combination of single 
finger device cells. Each single finger device cell is represented by a hierarchy of four 
models, which, in turn, are assembled together using models for interconnects to 
represent an arbitrary multifingered device cell, illustrated in Fig. 33. The four models 
are as follows: off mesa or boundary parasitic model; interelectrode parasitic model; 
on-mesa parasitic model and intrinsic model. 

[0176] The off-mesa parasitic model is illustrated in FIG. 34. This model 
represents the parasitics that exist outside the active FET region for each gate finger. 
In this model, the fringing capacitance of each gate finger off the active device region 
as well as the off-mesa gate finger resistance is modeled. 

[0177] The interelectrode parasitic model and corresponding equivalent circuit 
are illustrated in FIGS. 35-37. This model represents parasitics between the metal 
electrodes along each gate finger. The following fringing capacitance parasitics are 
modeled for the gate-to-source air bridge; drain-to-source air bridge; gate-to-source 
ohmic; gate-to-drain ohmic and source-to-drain ohmic as generally illustrated in FIG. 
36. 

[0178] The on-mesa parasitic model and corresponding equivalent circuit are 
illustrated in FIGS. 38 and 39. This model represents that parasitics around the active 
FET region along each gate finger including various capacitance fringing parasitics and 
resistive parasitics. In particular, the gate-to-source side recess; gate-drain- side recess; 
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gate-source access charge/doped cap; and gate-drain access charge/doped cap 
capacitance fringing parasitics are modeled. In addition, the gate metallization and 
ohmic contact resistive parasitics are modeled. 

[0179] The intrinsic model and corresponding equivalent circuit are illustrated 
in FIGS. 40 and 41. The intrinsic model represents the physics that predominately 
determine the FET performance. In particular, the DC and current voltage response 
can be determined by physics based analytical equations for magnitude and location 
of intrinsic charge which are generally know in the art, for example, as disclosed in 
"Nonlinear Charge Control in AlGaAs/GaAs Modulation-Doped FETs", by Hughes, 
et al., TEEE Trans. Electron Devices. Vol. ED-34, No. 8, August 1987. The small 
signal model performance is modeled by taking a derivative of the appropriate charge 
or current control equations to derive various terms such as RI, RJ, RDS, RGS, RGD, 
GM, TAU, CGS, CDS and CGD. Such control equations are generally known in the 
art and disclosed in detail in the Hughes, et al. reference mentioned above, hereby 
incorporated by reference. The noise performance may be modeled by current or 
voltage perturbation analysis "Noise Characteristics of Gallium Arsenide Filed-Effect 
Transistors" by H. Statz, et al. IEEE-Trans. Electronic Devices, vol. ED-21, No. 9, 
September 1974 and "Gate Noise in Field Effect Transistors at moderately High 
Frequencies" by A. Van Der Ziel, Pro. IEEE, vol 5 1 , March 1 963 "Gate Noise in Field 
Effect Transistors at Moderately High Frequencies" , by H. Statz, IEEE Trans. Electron 
Devices , vol. ED-21, No. 9, September 1974. "Noise Characteristics of Gallium 
Arsenole Field Effect Transistors", by Statz et al, TEEE Trans. Electron Devices, vol. 
ED-21, No. 9, September 1974. 

[01 80] An example of a parasitic model for use with the S-parameter microscopy 
discussed above is illustrated in FIGS. 42-49. Although a specific embodiment of a 
semiconductor device is illustrated and described, the principles of the present 
invention are applicable to various semiconductors devices. Referring to FIG. 42A, 
a Pi-FET is illustrated. As shown, the PiFET has four gate fingers. The four fingered 
Pi-FET is modeled in FIG. 42B. In particular, FIG 42B illustrates an equivalent circuit 
model for Pi-FET illustrated in FIG. 42A as implemented by a known CAD program, 
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for example, LIBRA 6.1 as manufactured by Agilent Technologies. As shown, the 
equivalent circuit models does not illustrate all of the equivalent circuit elements or 
network connections involved with implementing the parasitic embedding models, but 
rather demonstrates a finished product. The actual technical information regarding the 
construction of the network and its equivalent circuit elements are normally provided 
in schematic view. An important aspect of parasitic modeling relates to 

modeling of multi-gate fingered devices as single gate finger devices. As used herein, 
a single unit device cell refers to a device associated with a single gate finger. For 
example, a four fingered Pi-FET as illustrated in FIG. 42A is modeled as four unit 
device cells. 

[0181] Initially, the four finger Pi-FET illustrated in FIG. 42A, is modeled as a 
single finger unit device cell 100 with an intrinsic model 102, as shown in FIGS. 43 
and 44. In particular, the Pi-FET intrinsic FET model 104 is substituted for the block 
102 defining a first level of embedding. As shown in FIG. 44, the parameter values 
for the Pi-FET intrinsic model are added together with the parameter values for the 
single fingered unit device cell intrinsic model. The intrinsic device model 104 may 
be developed by S-parameter microscopy as discussed above. Next, as illustrated in 
FIG. 45, the interconnect layout parasitic elements are added to the equivalent model 
by simply adding the model terms to the value of the appropriate circuit element to 
form a single unit device cell defining a second level of embedding. Once the single 
unit device cell is formulated, this device is used to construct models for multi-fingered 
devices. In this case, a Pi-FET with four gate fingers is modeled as four single finger 
device unit cells as shown in FIG. 46. Subsequently, the off-mesa layout parasitic 
elements are connected to the multi-fingered layout, defining a third level of 
embedding as illustrated in FIG. 47. These off-mesa layout parasitic elements, 
generally identified with the reference numerals 1 08 and 1 1 0, are implemented as new 
circuit elements connected at key outer nodes of the equivalent circuit structure. 
Subsequently, a fourth level of embedding is implemented as generally illustrated in 
FIG. 48. In particular, an inductor model is connected to the sources of each of the 
various unit device cells to represent the metallic bridge interconnection, as generally 
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shown in FIG. 48. Lastly, as illustrated in FIG. 49, a fifth level of embedding is 
implemented in which the feed electrodes model 1 14 and 1 16 are modeled as lumped 
linear elements (i.e. capacitors inductors) as well as the distributive elements (i.e. 
microstrip lines and junctions) to form the gate feed and drain connections illustrated 
in FIG. 53. As shown, the distributive elements are distributed models for microstrip 
elements as implemented in LIBRA 6.1. 

EXTRACTION METHOD FOR UNIQUE DETERMINATION OF FET 
EQUIVALENT CIRCUIT MODELS 

[0182] The method for determining FET equivalent circuit parameters as 
discussed above is illustrated in FIGS. 50-55. This method is based on an equivalent 
circuit model, such as the common source FET equivalent circuit model illustrated in 
FIG. 5. Referring to FIG. 5 OA, a model is initially generated in step 122. In 
accordance with an important aspect of the algorithm, the equivalent circuit parameters 
are based upon measured FET S-parameters. Measurement of S-parameters of 
semiconductor devices is well known in the art. FIG. 53 A is a Smith chart illustrating 
exemplary measured S-parameters SI 1, S12 and S22 for frequencies between 0.05 to 
40 GHz. FIG. 53B represents a magnitude angle chart for the measured S-parameter 
S21 from frequencies from 0.05 to 40 GHz. After the S-parameters are measured, as 
set forth in step 124 (FIG. 5 OA), it is ascertained whether the measurements are 
suitable in step 126. This is either done by manually inspecting the test result for 
anomalies, or by algorithms to validate the test set. If the measurements are suitable, 
the S-parameter measurements are stored in step 128. 

[0183] A space of trial starting feedback impedance point values, for example, 
as illustrated in Table 9 is chosen. Then, a direct model extraction algorithm, known 
as the Minasian algorithm, is used to generate preliminary values for the equivalent 
circuit model parameters, for each value of starting feedback impedance. Such 
extraction algorithms are well known in the art, for example, as disclosed "Broadband 
Determination of the FET Small Equivalent Small Signal Circuit" by M. Berroth, et 
aL, IEEE - MTT , Vol. 38, No. 7, July 1990. Model parameter values are determined 
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for each of the starting impedance point values illustrated in Table 3. In particular, 
referring to FIG. 5 OA, each impedance point in Table 9 is processed by the blocks 130, 
132, etc. to develop model parameter values for each of the impedance point in order 
to develop an error metric, which, in turn, is used to develop a unique small signal 
device model, as will be discussed below. The processing in each of the blocks 130, 
132 is similar. Thus, only a single block 130 will be discussed for an exemplary 
impedance point illustrated in Table 9. In this example, the feedback impedance point 
17 which correlates to a source resistance ohm of L7Q and a source inductance L s 
of 0.0045pH is used. 



TABLE 9 

Trial Starting Feedback, Impedance Space Point Values 



Impedance 


Resistance 


Inductance 


Point 


(Rs) 


(Ls) 


1 


0.1 Q 


0.0045 pH 


2 


0.2 Q, 


0.0045 pH 


3 


0.3 Q 


0.0045 pH 


4 


0.4 Q 


0.0045 pH 


5 


0.5 Q 


0.0045 pH 


6 


0.6 Q 


0.0045 pH 


7 


0.7 a 


0.0045 pH 


8 


0.8 Q 


0.0045 pH 


9 


0.9 Q 


0.0045 pH 


10 


1.0 Q 


0.0045 pH 


11 


1.1 Q 


0.0045 pH 


12 


1.2 Q 


0.0045 pH 


13 


1.3 Q. 


0.0045 pH 


14 


1.4 a. 


0.0045 pH 


15 


1.5 Q 


0.0045 pH 


16 


1.6 Q 


0.0045 pH 


17 


1.7 Q 


0.0045 pH 


18 


1.8 Q 


0.0045 pH 


19 


1.9 a 


0.0045 pH 


20 


2.0 a 


0.0045 pH 


21 


2.1 Q 


0.0045 pH 


22 


2.2 Q 


I 0.0045 pH 
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23 


2.3 Q 


0.0045 pH 


24 


2.4 Q 


0.0045 pH 


25 


2.5 Q 


0.0045 pH 


26 


2.6 Q 


0.0045 pH 


27 


2.7 a 


0.0045 pH 


28 


2.8 a 


0.0045 pH 


29 


2.9 Q 


0.0045 pH 


30 


3.0 Q 


0.0045 pH 



[0184] For the selected value, R, = 1 .7 ohms, initial intrinsic equivalent circuit 
parameters and initial parasitic equivalent circuit parameter are determined, for 
example, by the Minasian algorithm discussed above and illustrated in Tables 10 and 
1 1 as set forth in steps 134 and 136. In step 138, the simulated circuit parameters are 
compared with the measured S-parameters, for example, as illustrated in FIGS. 54 A 
and 54B. Each of the processing blocks 130 and 132 etc. goes through a fixed number 
of complete cycles, in this example, six complete cycles. As such, the system 
determines in step 140 whether the six cycles are complete. 



TABLE 10 

Initial "Intrinsic" Equivalent Circuit Parameters 



Intrinsic Equivalent 


Initial Solution 


Circuit Parameter 




Cgs 


0.23595 pF 


Rgs 


91826 Q 


Cgd 


0.0177 pF 


Rgd 


100000 Q 


Cds 


0.04045 pF 


Rds 


142.66 Q 


Gm 


142.1025 mS 


Tau 


0.1 pS 
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TABLE 11 

Initial "Parasitic" Equivalent Circuit Parameters 



intrinsic Equivalent 


initial Solution 1 


Circuit Parameter 




Rg 


3.0 Q 


Lg 


0.014 nH 


Rs 


1.7 n 


Ls 


0.0045 nH 


Rd 


2.5 Q 


Ld 


0.024 nH 



[0185] Each cycle of the processing block 130 consists of a direct extraction 
followed by an optimization with a fixed number of optimization iterations, for 
example 60. By fixing the number of extraction-optimization cycles along with the 
number of optimization iterations, a fixed "distance" or calculation time which the 
model solution must be derived is defined. As such, the algorithm implements a 
convergence speed requirement of the global error metric by setting up an environment 
where each trial model solution competes against each other by achieving the lowest 
fitting error over a fixed calculation time thus causing "race" criteria to be 
implemented where "convergence speed" is implicitly calculated for each processing 
block 130, 132 etc. 

[0186] After the system determines whether the racing is done in step 140, the 
system proceeds to block 142 and optimizes model parameters. Various commercial 
software programs are available, for example, the commercially available, LIBRA 3.5 
software as manufactured by HP-eesof may be used both for circuit simulation as well 
as optimizing functions. The optimization is performed in accordance with the 
restrictions set forth in Table 12 with the addition of fixing the feedback resistance R s 
to a fixed value. 
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TABLE 12 

Environment Used for Competitive Solution Strategy, as Implemented in this 

Example 



Implementation Parameter 




Circuit Simulator and Optimizer 


Libra 3.5 


Optimization Algorithm 


Gradient 


Optimization Error Metric 


Mag and angle of 
S11,S21,S12,andS22 
from 4 to 40 GHz 


Number of Iterations 


60 


Number of Extraction/Optimization 
Cycles 


6 



[0187] By fixing the value for R, this segment of the algorithm confined to 
creating a trial model solution for only the trial feedback impendence point with which 
it started. Table 13 illustrates the optimized intrinsic equivalent parameter values 
using commercially available software, such as LIBRA 3.5. These values along with 
the optimized parasitic values, illustrated in Table 14, form the first optimized model 
solution for the first extraction-optimization cycle (i.e. one of six). The optimized 
model parameters are then fed back to the function block 134 and 136 (FIG. 50A) and 
used for a new initial model solution. These values are compared with the measured 
S-parameter value as illustrated in FIGS. 54A and 54B. The system repeats this cycle 
for six cycles in a similar fashion as discussed above. After the six extraction- 
optimization cycle, the final trial model solution for the trial impendence point 17 is 
complete along with its final fitting error to the measured data to form the new error 
metric 144. In accordance with an important aspect, the extraction-optimization 
algorithm makes the final optimization fitting error for each point implicitly carry 
information about both the measured to model fitting error and the speed of 
convergence. It does so by the fixed optimization time constraint which sets up a 
competitive race between the various trial model solutions. 
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TABLE 13 

Optimized "Intrinsic" Equivalent Circuit Parameters 



Intrinsic Equivalent 


initial Solution 


Circuit Parameter 




Cgs 


0.227785 pF 


Rgs 


65247 Q. 


Cgd 


0.017016 pF 


Rgd 


1 30820 Q 


Cds 


0.047521 pF 


Rds 


160.18 Q 


Gm 


135.74 mS 


Tail 


0.446 pS 



TABLE 14 

Optimized "Parasitic" Equivalent Circuit Parameters 



Intrinsic Equivalent 
Circuit Parameter 


Initial Solution 


Rg 


4.715 Q 


Lg 


0.02903 nH 


Rs* 


1.7 n 


Ls 


0.002102 nH 


Rd 


3.2893 Q 


Ld 


0.0317 nH 



[01 88] The implementation of the extraction optimization cycles makes the best 
and fastest solving solution appear as a global minima for the final fitting error in step 
146 of all of the trial impedance points as generally shown in FIGS. 5 1 and 52. More 
specifically, referring to FIG. 5 1 the global minima solution using the new error metric 
is found around R s =l .7 ohms. Tables 15 and 16 list the final model equivalent circuit 
parameters for this global solution, including the intrinsic and parasitic parameter as 
set forth in step 148 (FIG. 50B). 
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TABLE 15 

Global Solution for "Intrinsic" Equivalent Circuit Parameters 



1 Intrinsic Equivalent Circuit 


Initial Solution 


Parameter 




Cgs 


0.227745 pF 


Rgs 


64242 Q 


Cgd 


0.017019 pF 


Rgd 


133450 Q 


Cds 


0.047544 pF 


Rds 


160.1791 Q 


Gm 


135.7568 mS 


Tau 


0.443867 pS 



TABLE 16 

Global Solution "Parasitic" Equivalent Circuit Parameters 



Extrinsic Equivalent Circuit 


Initial Solution 


Parameter 




Rg 


4.711895 Q 


Lg 


0.029314 nH 


Rs 


1.7 n 


Ls 


0.002104 nH 


Rd 


3.309899 Q 


Ld 


0.031671 nH 



[0189] In order to test the accuracy of the solution, the final model for solutions 
are compared with the measured S-parameter values as shown in FIGS. 55 A and 55B. 
As shown, there is good correlation between the simulated model values and the 
measured S-parameters values thus verifying that the simulated model values represent 
a relatively accurate and unique small signal device model. 

[0190] Obviously, many modifications and variations of the present invention are 
possible in light of the above teachings. Thus, it is to be understood that, within the 
scope of the appended claims, the invention may be practiced otherwise than as 
specifically described above. 
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[0191] What is claimed and desired to be covered by a Letters Patent 
follows: 
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